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1.  SUMMARY 


Utilization  of  the  full,  seismic-waveform  based  techniques  for  discrimination  purposes  is 
currently  effective  and  efficient  for  large  magnitude,  >  M  4.5,  well-recorded  events.  However, 
a  need  to  reduce  the  magnitude  threshold  requires  predictive  capabilities  at  shorter  periods, 

<~10  seconds,  if  waveform  based  techniques  are  to  be  used  in  routine  event  evaluations. 

Recent  efforts  to  improve  the  seismic  waveform,  prediction  capabilities  of  existing  models  have 
proven  successful  by  reducing  an  acceptable  period  for  prediction  to  15-20  seconds.  As 
models  begin  to  accurately  predict  shorter  periods,  they  will  begin  to  incorporate  higher  order 
structures,  i.e.  larger  and  sharper  wave  speed  contrasts,  and  they  will  also  be  effective  in  source 
discrimination  at  shorter  periods  and  high  noise  environments.  However,  as  the  minimum 
period  continues  to  drop,  effects  of  attenuation  and  anisotropy  will  become  more  significant  and 
require  quantification.  As  such,  current  generation  models  lack  both  heterogeneous  anelastic 
attenuation  and  a  more  robust  treatment  of  anisotropy.  We  plan  on  pushing  the  predictive 
capabilities  of  the  current  3D  Middle  East  wave  speed  model  to  shorter  periods.  In  doing  so, 

3D  anleastic  attenuation  and  anisotropic  models  will  be  incorporated  into  a  seamless  seismic 
Earth  model  of  the  Middle  East  region.  Anelastic  attenuation  will  be  embedded  along  side  a  3D 
wave  speed  model  and  the  anisotropic  component  will  be  parameterized  as  transverse  isotropy 
with  a  variably  oriented  symmetry  axis. 

Achieving  these  ambitious  goals  requires  a  waveform  based,  tomographic  inversion  capable  of 
handling  seismic  waveform  propagation,  3D  variations  in  attenuation,  and  general  anisotropy. 
As  with  previous  efforts  use  of  the  SPECEEM3D  Globe  and  its  adjoint  tomographic  capabilities 
will  perform  the  necessary  forward  and  time-reversed  calculations.  Demands  required  by  the 
short  period  requires  use  of  large  clustered  computer  systems  and  improved  /  state-of-the-art 
computational  units  to  efficiently  simulate  the  seismic  wavefield.  An  initial  step  of  source 
characterization  demonstrated  that  this  is  invaluable  to  accurately  predict  waveform  amplitudes 
with  minimal  change  to  the  source  moment  tensor,  but  only  at  stations  near  amplitude  nodes. 
Predicting  wave-field  amplitudes  are  increasingly  important,  as  anelastic  attenuation  is  highly 
dependent  on  variations  in  seismic  amplitudes.  Initial  inversions  for  seismic  wave  speeds  began 
at  longer  periods  and  later  inversions  incorporated  additional  shorter  period  energy.  Our  results 
show  wave  speed  anomalies  began  to  localize  and  strengthen  in  amplitude;  similar  effects  for 
attenuation  and  anisotropy  are  also  expected.  The  final  model,  with  a  enhanced  predictive 
capabilities,  3D  anelastic  attenuation,  and  anisotropy  with  a  variably  oriented  symmetry  axis 
will  be  a  major  advance  for  the  Nuclear  Explosion  Monitoring  community  for  event 
discriminate  at  smaller  magnitudes. 

2.  INTRODUCTION 

This  work  aims  to  improve  predictions  of  full-waveform  seismic  signals  from  regional  events 
within  the  Middle  East  region.  When  complete,  responses  from  events  within  the  study  region 
will  be  characterized  and  predicted  at  shorter  periods.  Previous  results  utilizing  only  wave 
speeds,  Vp,  Vsv  and  Vsh,  were  able  to  adequately  reproduce  recorded  events  down  to  15 
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seconds.  Updates  to  the  3D  wave  speed  model  significantly  reduced  travel  time  misfit  but  did 
not  improve  differences  in  amplitude.  Further  enhancements  will  require  the  incorporation  of  3D 
attenuation  and  more  complex  treatment  of  anisotropy.  Iterative  improvements  to  the  3D  Earth 
model,  including  wave  speed,  attenuation  and  anisotropy,  are  accomplished  through  an  adjoint 
method  that  identifies  regions  of  the  model  that  require  improvement.  The  adjoint  method 
incorporates  the  physics  of  wave  propagation  and  finite-frequency  effects  and  can  improve 
predictions  of  travel  times  and  amplitudes,  i.e.  the  full  waveform. 

Improvement  of  the  3D  Earth  model  will  radically  enhance  source  discrimination 
investigations.  Significant  limitations  to  full-waveform  source  studies  are  the  lack  of  Earth 
models  that  accurately  predict  the  path-specific  Green’s  functions.  While  these  source 
investigations  are  quite  robust  at  discriminating  between  a  variety  of  source  types  (e.g.  Eiu  et  al, 
2004;  Minson  and  Dreger,  2008;  Eord  et  ah,  2009;  Kim  et  ah,  2011),  the  limited  predictive 
capabilities  of  wave  speed  models  restrict  their  periods  use  to  greater  than  30  -  50  seconds.  Use 
of  these  long  periods  necessitates  sources  to  be  a  minimum  magnitude,  nominally  >  4.5  Mw. 
Efforts  to  enhance  the  waveform  predictive  capabilities  of  wave  speed  models  have  recently 
begun  to  shift  the  minimum  period  to  --15  seconds  at  regional  distances.  The  recent  and  future 
improvement  of  these  full  waveform  predictive  models  will  allow  the  investigation  of  sources  at 
even  lower  magnitudes,  ~  3.0  Mw,  dependent  on  epicentral  distance. 


10"  20'  30"  40"  50"  60"  70"  80"  90"  100" 


Figure  1.  Map  of  the  Middle  East  Region.  Map  of  the  study  region  including  events 
as  red  circles  and  stations  as  blue  triangles.  The  original  mesh,  black  outline,  was 
approximately  45  degrees  on  a  side  and  easily  capable  of  computing  synthetics  to  a 
minimum  period  of  11  seconds.  A  new  expanded  mesh  is  90  degrees  on  a  side  and 
includes  the  many  more  events  and  stations  not  included  in  the  original  mesh  design. 
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3.  TECHNICAL  APPROACH 


3.1.  Earth  Model  Tomography 

Extending  the  predietive  period  range  of  the  Middle  East  model  will  be  undertaken  using  the 
adjoint  eapabilities  of  SPECEEM3D_GEOBE.  Essential  to  the  adjoint  tomographie  method 
(Tromp  et  ah,  2005)  and  improving  3D  Earth  models  is  the  ability  to  calculate  the  regular  and 
adjoint  wavefields  in  an  existing  3D  Earth  model.  The  regular  wavefield,  or  forward  field,  is 
defined  as  the  displacement  at  all  points  for  all  times  from  a  traditional  seismic  source.  While  the 
adjoint  wavefield  is  defined  as  the  displacement  at  all  points  for  all  reversed  times  from  adjoint 
waveform  sources.  The  adjoint  waveform  source  is  a  user-defined  measure  between  the  true 
response  at  a  seismic  station,  the  data,  and  the  predicted  response,  the  synthetic.  This  measure 
can  be  a  simple  difference  between  the  data  and  synthetics,  a  least  squares  differences,  a  metric 
isolating  travel  time  or  amplitude  differences,  or  any  combination.  Eor  this  specific  project,  we 
use  a  multi-taper  methodology  to  measure  frequency-dependent  delays  in  travel  time  and 
differences  in  amplitude.  A  summation  of  the  regular  and  adjoint  fields  over  time  generates 
Erechet  derivatives  that  isolate  deficiencies  in  the  current  model,  i.e.,  where  the  model  requires 
improvement.  The  adjoint  method  implemented  here  sums  the  regular  and  adjoint  fields  during 
back-propagation,  or  while  the  forward  field  runs  backwards  in  time.  This  summation  eliminates 
the  requirement  of  storing  the  regular  and  adjoint  fields  for  all  times. 

Adjoint  techniques  use  full  seismic  wave  propagation  physics  and  incorporate  finite- 
frequency  wave  propagation  effects.  These  methods  have  been  used  with  impressive  results  in 
the  Middle  East  (Savage  et  ah,  2011,  2012,  2014),  Southern  California  (Tape  et  ah,  2009),  and 
Europe  (Zhu  et  ah,  2013).  Many  of  the  current  advanced  seismic  tomographic  techniques 
include  various  fundamentals  of  the  adjoint  methodology,  but  do  not  exploit  the  full  waveform  to 
improve  the  predictive  response  of  surface  and  body  waves  in  tandem.  A  simple  example  of 
creating  the  kernels  is  displayed  in  Eigure  2  with  a  single  source  and  single  receiver. 

This  new  adjoint  tomography  experiment  will  increase  the  seismic  data  by  '-'30%  through  the 
inclusion  of  more  recent  seismic  events  (Eigure  1).  The  data  set  will  also  vastly  increase  due  to 
the  ability  to  model  smaller  magnitude  events  at  shorter  periods.  Event  and  station  distributions 
(Eigure  1)  have  not  changed  much  from  the  initial  experiment,  but  regions  of  interest  tectonically 
have  instrumentally  expanded  thus  adding  more  seismic  data. 

Efforts  to  extend  the  period  range  to  even  shorter  periods  require  higher  resolution  images  of 
wave  speeds,  plus  incorporation  of  increasingly  important  effects  of  anelastic  attenuation,  Q'^ 
and  anisotropy.  With  decreasing  period  these  second-order  effects,  at  long  period  <30  seconds, 
become  important  and  can  radically  alter  seismic  paths  and  recorded  responses.  The  goal  of  this 
effort  is  to  further  improve  models  of  seismic  wave  speeds.  Vs  and  Vp,  to  predict  seismic 
waveforms  at  periods  approaching  2  seconds.  In  addition,  anelastic  attenuation  and  anisotropy 
will  also  be  accounted  for  through  improved  use  of  amplitude  anomalies  and  parameterization. 
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Figure  2,  Example  of  adjoint  kernel  creation.  Example  of  the  adjoint  tomography 
method,  where  time  increases  in  each  row  starting  from  the  bottom,  the  first  two 
columns  are  the  regular  and  adjoint  wavefields,  the  third  column  is  the  interaction, 
or  multiplication,  between  the  first  two  columns  and  the  fourth  column  is  the  kernel, 
or  time  summation/integration  of  the  third  column.  Notice  the  regular  wavefield 
propagates  from  the  source,  star,  to  the  receiver,  box,  in  forward  time  (bottom  row 
to  top  row)  and  the  adjoint  field  propagates  from  the  receiver  to  the  station,  but  in 
reverse  time  (top  row  to  bottom  row).  From  Tromp  et  al.  (2005) 


3.2.  Source  Re-evaluations 

In  addition  to  improvement  to  the  Earth  model,  seismic  sources  are  reevaluated  in  an  effort  to 
reduce  artifacts.  Source  related  artifacts  in  the  final  model  are  a  result  of  errors  in  the  fault  plane 
solution  and  hypocenter.  We  use  a  straightforward  moment  tensor  inversion  outlined  in  Eiu  et  al. 
(2004)  to  account  for  and  eliminate  these  errors.  The  advantage  of  the  methodology  is  the 
flexibility  to  use  synthetic  seismograms  from  any  methodology,  including  the  improved  3D 
Earth  model.  We  solve  for  the  best  fitting  moment  tensor  solution  and  depth,  constrained  by  a 
zero-trace  solution.  Moreover,  we  also  solve  for  the  best  fitting  double  couple,  full  moment 
tensor,  and  perform  an  error  analysis.  In  practice,  all  moment  tensor  solutions  are  quite  similar 
to  the  preferred  solution  and  to  the  initial  solution  provided  by  either  the  Global  CMT  solution  or 
previous  work  by  Covellone  and  Savage  (2012).  Solution  errors  are  normally  quite  small  as 
defined  by  the  spread  of  P-  and  T-axes,  and  improvements  to  the  sources  show  up  at  azimuths 
close  to  nodal  planes  of  the  radiation  pattern.  Recorded  responses  show  rapid  changes  in 
amplitudes  around  nodal  planes,  thus  small  changes  in  the  source  solution  can  dramatically  alter 
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a  the  response  for  a  near  nodal  station.  With  these  minor  changes,  we  anticipate  a  further 
increase  in  the  double  couple  component  for  all  events  as  the  model  continues  to  improve. 

The  inversion  requires  synthetics  from  initial  solution,  obtained  from  Global  CMT  or 
Covellone  and  Savage  (2012),  a  perturbed  depth  source,  and  sources  for  each  of  the  six  moment 
tensor  elements.  Each  of  these  sources  is  computed  using  the  current  iteration  model,  M14,  thus 
requiring  8  simulations  per  source  and  a  short  hiatus  of  the  tomography  iterations.  However, 
once  the  simulations  are  complete,  source  inversions  are  very  fast  and  exploited  during  a 
bootstrap,  error  analysis. 


Vertical 
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_ teration _ 

Figure  3.  Misfit  between  data  and  synthetics  against  iteration.  Total  misfit 
(bottom)  is  separated  into  individual  components  (columns  in  first  and  second  rows) 
and  into  short  (top  row)  and  long  period  (middle  row)  bands.  The  short  period  band 
involves  primarily  body  waves,  while  the  long  period  band  is  for  surface  waves.  A 
step  in  the  misfit  at  iteration  4  is  due  to  a  change  in  both  period  bands. 


4.  RESULTS  AND  DISCUSSION 

Adjoint  tomographic  inversions  were  completed  up  to  iteration  14  and  include  transverse 
isotropy  with  a  vertical  symmetry  axis  using  measurements  of  phase  delay.  Using  the  flexible 
measurement  tool  FLEXWfN  (Maggi  et  al.  2009),  measurements  of  misfit  were  easily  made  for 
surface  and  body  waves  simultaneously  on  three  components  and  in  two  different  period  bands. 
The  misfit  for  all  components  and  the  total  misfit  are  displayed  in  Figure  3. 
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During  each  iteration,  improvements  to  the  model  are  made  and  the  misfit  is  reduced  for  all 
components  and  all  period  bands.  The  most  significant  improvements  are  during  the  initial 
iterations,  but  as  the  iterations  progresses  the  misfit  and  the  improvements  to  the  model  are 
reduced  and  the  misfit  plateaus,  both  are  expected  behaviors.  A  significant  step  in  the  misfit  at 
iteration  4  is  due  to  a  modification  of  the  period  bands  for  both  long  and  short  periods.  The  long 
period  was  expanded  from  a  range  of  40-120  seconds  to  30-120  seconds,  but  the  short  period 
was  reduced  from  14-40  seconds  to  14-30  seconds.  The  reduction  in  the  short  period  band 
accounts  for  most  of  the  drop  in  misfit  at  this  iteration.  A  second  change  to  the  parameterization 
was  made  at  iteration  7  where  the  lateral  smoothing  was  reduced  from  100  km  to  50  km,  but  this 
appears  to  have  little  impact  on  the  overall  improvement  of  the  model. 

At  iteration  14,  the  period  ranges  are  currently  14-30  seconds  and  20-120  seconds.  Iteration 
14  also  represents  the  transition  to  a  larger  computational  and  model  domain  that  includes  many 
more  events,  stations  and  usable  data  than  in  previous  iteration.  Figure  4.  In  the  new  model,  for 
regions  outside  of  the  original  model  domain  the  initial  wave  speed  model,  MOO,  was  used  while 
model  iteration  13,  Ml  3,  was  used  within.  This  expansion  of  the  model  does  not  significantly 
alter  the  total  misfit  between  M13  and  M14.  The  increase  in  events,  stations,  and  ray-based  path 
coverage  is  displayed  in  Figure  4.  Expanding  the  model  domain  introduces  stations  in  northern 
Asia,  eastern  Europe,  and  in  the  Indian  Ocean  along  with  a  vast  number  of  events  in  Tibet  and 
post-2007  events.  The  path  coverage  represents  usable  data  with  measurement  windows  as 
determined  by  ELEXWIN. 


10'  20'  30'  40'  50'  60'  70'  80'  90'  100'  10'  20'  30’  40'  50’  60’  70’  80'  90’  100’ 

Figure  4,  Comparison  of  path  coverage  between  old  and  new  models.  Ray  path 
coverage  in  the  original  model  domain  (left)  and  the  new,  larger  model  domain 
(right).  The  original  model  domain  contained  182  events  and  356  stations,  but  not  all 
the  stations  had  workable  signals.  The  new  model  domain  includes  468  events  and 
241  stations,  all  with  usable  signals. 


The  current  model  for  the  Middle  East,  Ml 4,  is  displayed  in  Eigure  5  as  absolute  wave  speed 
and  a  wave  speed  perturbation.  The  wave  speed  perturbation  indicates  a  total  change  of  around 
6%  in  vertical  shear  wave  speed,  Vsv  from  the  initial  3D  model,  MOO,  S2.9EA  (Kustowski  et  ah, 
2008).  A  number  of  well  known,  large-scale  tectonic  features  appear  within  the  model  further 
suggesting  that  the  inversion  is  converging  on  a  true  representation  of  the  Earth.  Beneath  the 
Persian  Gulf,  there  is  a  subducted  Neo-Thethys  oceanic  lithosphere,  faster  than  average,  trending 
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to  the  northeast  and  extending  throughout  mueh  of  the  upper  mantle.  Low  wave  speeds  in  the 
eentral  and  southern  Red  Sea  extending  into  the  Gulf  of  Aden  are  quite  signifieant  down  to  at 
least  200  km  depth.  Finally,  subdueting  Indian  lithosphere  is  identified  beneath  the  western  most 
edge  of  Tibet  as  well  as  thiekened  crust  not  accounted  for  in  the  initial  model  appears  as  a 
reduction  in  wave  speed  at  shallow  depths,  <  75km. 


dvswsv 

Figure  5,  Current  iteration  wave  speed  model  of  the  Middle  East.  Absolute 
vertical  shear  wave  speed,  Vsv,  is  shown  in  the  left  column  and  wave  speed 
perturbation  of  Vsv  is  in  shown  the  right  column  plotted  at  depths  of  75,  100,  150 
and  200  km. 
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Figure  6,  Phase  and  amplitude  misfits  for  body  waves.  Histograms  of  phase,  dT, 
and  amplitude,  dlnA,  mismatch  from  the  initial  model,  MOO,  and  the  current  model. 
The  initial  model  is  the  single  line  in  red  and  the  current  model  is  the  filled  light  blue 
bars.  Improvement  to  both  the  phase  and  amplitude  are  apparent  across  all 
components.  Phases  are  body  waves  measured  between  14  and  30  seconds  period. 


Improvement  in  the  predictive  capabilities  of  the  current  iteration  model  is  demonstrated  in 
Figures  6  and  7.  For  the  short  period  band,  primarily  body  waves,  the  phase  and  amplitude 
mismatch  is  shown  for  the  initial  model,  MOO,  and  the  current  model,  M14.  For  a  perfect  fit, 
phase  and  amplitude  mismatch  should  be  0.0.  Phase  measurements  for  ail  three  components 
show  a  vast  improvement  from  the  initial  to  the  current  model.  The  phase  misfit  distribution  is 
narrower  with  a  larger  peak  around  zero,  indicating  a  larger  portion  of  body  wave  arrivals  are 
matched  with  zero-lag  time.  Results  are  similar  for  the  longer  period  band  that  highlights  surface 
wave  arrivals. 

Amplitude  misfits  are  also  reduced  while  only  measurements  of  phase  delay  are  used  to 
improve  the  model.  Improvements  to  amplitude  data  are  not  as  significant  as  those  for  phase 
delay,  but  this  is  expected.  The  combination  of  the  minimal  improvement  to  the  amplitude 
measurements  and  the  misfit  at  later  iterations  suggests  that  amplitude  measurements  should  be 
included  in  the  inversion  in  an  effort  to  update  the  attenuation  structure.  This  initially  will 
increase  the  misfit  due  to  the  added  measurements,  but  we  should  obtain  a  robust,  combined 
wave  speed  and  attenuation  model.  We  plan  on  starting  iterations  including  attenuation  during 
the  initial  stages  of  the  upcoming  year. 

An  example  of  improvements  to  the  model  and  its  predictive  capabilities  for  an  individual 
event-station  path  is  shown  in  Figure  8.  A  synthetic  response  on  the  western  side  of  the  model 
from  an  event  near  the  western  syntaxis  of  Tibet  demonstrates  improvements  to  the  model.  For 
ail  three  components  there  is  a  significant  improvement  in  the  fit  to  the  data.  Surface  waves, 
both  Rayleigh  and  Love,  show  good  travel  time  agreement  using  the  current  model  as  opposed  to 
large  time  shifts  seen  from  responses  from  the  initial  model.  Higher  frequency  body  waves. 
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particularly  on  the  radial  component  at  712  seconds,  show  improved  agreement  with  the  data  in 
the  updated  model.  These  radial  S  wave  arrivals  appear  quite  poorly  matched  when  using  the 
initial  model,  in  amplitude  and  phase.  Response  of  the  current  model  for  these  arrivals  shows  a 
dramatic  improvement  in  phase  but  at  the  expense  of  increased  amplitude.  This  might  indicate  a 
need  for  greater  attenuation  along  these  paths.  As  shown  by  the  histograms  in  Figures  6  and  7, 
many  of  the  synthetic  responses  show  improvement  in  a  similar  manner  as  demonstrated  for  the 
single  event-station  path  in  Figure  8. 
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Figure  7,  Phase  and  amplitude  misfits  for  surface  waves.  Same  as  Figure  8  except 
for  surface  waves  measured  between  24  and  120  seconds  period. 


Due  to  the  reduced  misfit  vs.  iteration  for  the  last  few  iterations  and  the  less  than  optimal  fit 
of  waveform  amplitudes,  demonstrated  by  the  raw  waveforms  and  the  amplitude  misfit  for  the 
entire  data  set,  upcoming  iterations  of  the  Middle  East  model  will  use  the  expanded  model  shown 
in  Figures  1  and  4.  This  modification  of  the  inversion  parameters  will  reduce  the  overall  misfit, 
but  is  required  as  further  iterations  using  the  current  parameters  and  data  sets  are  not  able  to 
account  for  the  variance  identified  in  the  data  (see  the  amplitude  mismatch  in  the  waveforms 
displayed  in  Figure  8).  New  major  changes  will  be  1)  an  introduction  of  amplitude 
measurements  and  the  ability  to  solve  for  a  3D  attenuation  model,  2)  the  relaxing  of  constraints 
on  anisotropy  by  solving  for  the  symmetry  axis  orientation  along  with  Vsv  and  Vsh,  and  3)  the 
introduction  of  a  vastly  larger  data  set  and  model  domain.  The  addition  of  a  3D  attenuation 
model  and  amplitude  measurements  have  successfully  been  implemented  in  a  European  model, 
building  on  prior  work  by  Zhu  et  al.  (2013).  Relaxing  the  anisotropy  constraints  will  improve 
fits  to  surface  wave  arrivals  for  continental  paths  with  complex  tectonic  histories. 
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5.  CONCLUSIONS 


Initial  work  improving  the  seismic  waveform  predictive  capabilities  for  the  Middle  East 
adjoint  tomography  model  demonstrated  that  enhancements  were  possible  up  to  15  seconds. 
These  results  at  longer  period  are  in  agreement  with  previous  studies  of  the  region.  Full 
waveforms  and  adjoint  techniques  were  used  to  identify  deficiencies  in  the  wave  speed  model. 
Initially,  to  progress  to  a  minimum  period  of  15  seconds,  transverse  isotropy  with  a  vertical 
symmetry  axis  was  used  in  the  mantle  to  stabilize  the  inversion.  Travel  times  of  all  types  of 
seismic  waves  showed  the  largest  variance  reduction,  whereas  there  was  minimal  improvement 
of  amplitude  variations.  Moreover,  use  of  3D  wave  speed  models  in  source  studies  showed  a 
dramatic  improvement  in  the  percent  double  couple  of  moment  tensors  over  the  entire  data  set. 

Future  work  will  continue  to  improve  the  seismic  waveform  predictive  capability  of  the 
Middle  East  adjoint  tomography  model.  We  intend  to  reduce  the  minimum  period  as  much  as 
possible  allowing  for  use  of  shorter  period  arrivals,  smaller  sized  sources,  and  data  with  smaller 
signal  to  noise  ratios.  To  reach  shorter  periods,  the  3D  anelastic  attenuation  variations  will 
become  increasingly  important  and  will  be  incorporated  into  the  Middle  East  model.  At  shorter 
periods,  knowledge  of  the  3D  anelastic  attenuation  structure  will  improve  the  fit  to  the 
amplitudes.  Incorporation  of  3D  anisotropy  will  also  be  essential  as  the  minimum  period 
continues  to  decrease.  We  will  relax  the  current  constraint  of  a  vertical  symmetry  axis  for 
transverse  isotropy  in  an  attempt  to  recover  the  anisotropy  in  the  Asian  continental  lithosphere 
and  surrounding  microplates. 

Improved  methods  for  routinely  computing  seismic  waveforms  at  shorter  periods  will  require 
larger  computational  and  clustered  systems,  and  improved,  highly  optimized  algorithms  and 
routines  to  reduce  unnecessary  computation. 
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Figure  8,  Comparison  of  waveforms  fits  for  a  single  event-station  path.  Single 
station-event  response  difference  between  the  initial  model,  left,  and  the  current 
model,  right.  Data  is  shown  in  black  and  synthetic  responses  are  in  red.  All  three 
components,  vertical,  radial  and  tangential  are  display  at  20  to  120  seconds.  The 
improvement  to  the  model  aligns  the  synthetic  response  with  the  data,  primarily  the 
surface  waves. 
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